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Abstract 

The probability distribution function for thermodynamics and econo- 
physics is obtained by solving an equilibrium equation. This approach 
is different from the common one of optimizing the entropy of the sys- 
tem or obtaining the state of maximum probability, which usually 
obtains as a result the Boltzmann distribution. The Gamma distri- 
bution is proposed as a better equation to describe the blackbody 
radiation in substitution of Planck's radiation equation. Also, a new 
form of entropy is proposed, that maintains the correct relation with 
the Clausius' formula. 
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1. INTRODUCTION 

Kirchhoff introduced the concept of a perfectly black body, being one 
which absorbs all incident radiation. He was the first in presenting the enunci- 
ation of the principle of the balance between emission and absorption around 
1860. The principle of equilibrium of emission and absorption between differ- 
ent agents is, essentially, the method that will be used in the present paper 
to develop a new equation that describes the radiation distribution of a black 
body. 

By 1900 Planck devised the known equation that describes the distribu- 
tion of radiation emitted by a black body, based on hunch and pure guessing. 
Although it is rather accurate over the entire range of frequencies it is sug- 
gested here that it is not appropriate to describe the blackbody radiation 
distribution. 
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"Although the introduction of energy quanta a century ago has led to the 
currently accepted explanation within quantum theory, there is still no firm 
conclusion as to whether or not blackbody radiation can be explained within 
classical physics." [7j 

One of the basic assumptions in this paper is that energy exists in packets, 
or quanta, which have continuous values. 

This deserves some explanation because this is against one of the most 
basic and accepted tenets of Physics, in the so-called area of quantum theory. 

The elementary particles called fermions, such as the electron and proton, 
have some very specific amounts of energy and also the orbits of the electrons 
in the atoms seem to be at discrete levels. Consequently, it can be safely 
accepted that those elementary particles satisfy the classical quantum theory. 
Moreover, it can be accepted that, when the electrons jump between levels 
in the atomic orbits, they emit and absorb photons with discrete levels of 
energy. 

However, this discrete character is not logically and necessarily extrapo- 
lated to the photons produced, for example, in the blackbody or an antenna 
radiation. It has not been proved that such photons have discrete levels of 
energy. On the contrary, electromagnetic radiation seems to have a continu- 
ous spectrum in the whole range of frequencies or wavelengths. The present 
author has not been able to find any experiment supporting the opposite 
position. 

It seems that some physicists might have some misconceptions about 
the current blackbody radiation theory. For example, the peer reviewers of 
Physica A, under the orders of Prof. C. Tsallis, used the argument that 
"It is obvious that the energy spectrum given by the Planck distribution is 
continuous", to reject a former version of the present paper. 

On the light of the points of view suggested by the following reference, 
such argument gives the impression of being wrong. 

Let us refer as "A&F" to the book of Alonso & Finn, ^^Fundamental Uni- 
versity Physics, Volume III, Quantum and Statistical Mechanics" , Addison- 
Wesley, 1968, which is referred to in the present paper as 

By definition. Quantum Mechanics deals with energy levels. This is sug- 
gested in Fig. 2-10 of A&F. 

In section 10-3, A&F carry out the proof of the "Maxwell-Boltzmann 
distribution." In particular, in the second paragraph of that section, the 
authors mention "i/ie number of particles present in each of the energy levels" . 
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Fig. 10-1, appearing in the same section, is also a good illustration of this 
concept. The derivation of the "probability distribution according to the 
Maxwell-Boltzmann distribution", equation (10.8), is very clear and should 
be reviewed to confirm that the authors A&F always use discrete energy 
levels. After using a simple mathematical optimization technique they find 
the "partition of maximum probability" and the law of "Maxwell-Boltzmann 
distribution", equations (10.9) and (10.11), which provide the integer number 
of particles in each energy level. 

It is obvious that the products of integer numbers by discrete energy levels 
produce discrete levels for the energy distribution. 

More relevant to the present paper is the derivation of the Planck radia- 
tion energy distribution, eq. (13.20), which is carried out by A&F in section 
13.6, that they call "maybe the most important application of the Bose- 
Einstein statistics" (rather similar to the Maxwell-Boltzmann distribution, 
given the additional considerations identified in the cited book). The authors 
mention that it is only an approximation that the energy spectrum could be 
considered continuous, under the condition that "the cavity be great relative 
to the mean radiation wavelength because the energy difference among two 
successive energy levels is sraall." This last condition is satisfied because, in 
general, we have a "big box", as pictured in Fig. 2-10 (b), compared with the 
small values obtained for energy, which is just the product of the frequencies 
of the packets by the Planck constant. But this does not mean that Quantum 
Mechanics does not consider energy levels or that those levels have somehow 
disappeared during the derivation. 

Moreover, the whole mentioned derivation produces the mathematical 
forms of the Planck radiation distributions, which are different from the dis- 
tribution obtained under the assumption of continuous energy values for the 
photons, as is proved in this paper. 

Some investigators have recognized the fact that there are two different 
maxima for the Planck function, which they justify using explanations that 
the present author and even their students consider implausible. Read, for 
example, the anecdotes after equations (1.24) and (1.25) of section 1.2.1 in 
the book by Bohren & Clothiaux [S]. Those authors conclude that "Failure 
to recognize that the maximum of a distribution function depends on how it 
is plotted has led and no doubt will continue to lead to errors." On the other 
hand, Soffer & Lynch [221 direct their quest to the spectral sensitivity of the 
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human eye in order to unify (within the eye) the different maxima, with the 
premise that the problem "results from the choice of units in which the solar 
spectrum is plotted." 

The same or a very similar argument seems to be getting out of con- 
trol with other investigators [IHl^Eni that write a variety of equations that 
supposedly make an equivalent "bookkeeping" of the "Wien peak" . For ex- 
ample ^3] reports at least five equations and refers to four of them saying 
that "[a] 11 four of these dispersion rules, with their corresponding graphs and 
formulas, represent the same physical spectrum." The same author, making 
reference to a publication of the American Association of Physics Teachers 
that reports quibbles, misunderstandings and mistakes in Physics, says that 
"The analysis and interpretation of continuous spectra has a reputation for 
leading people astray". However, the physicists and the Physics teachers in 
particular would not move a finger to correct those mistakes. The attempt 
to publish a former version of the present paper received the following reply, 
supposedly stated in the Statement of Editorial Policy: "AJP [a publication 
of the American Association of Physics Teachers] may only consider for publi- 
cation manuscripts dealing with concepts accepted by the general community 
of physicists" . 

There is no doubt that the current concepts about blackbody radiation, 
Planck equations, Wien displacement law, Boltzmann entropy and others 
that are challenged in the present paper are unanimously accepted by the 
scientific community. Consequently, the absolute acceptance of those (and 
other) concepts puts an unsurmountable barrier to the possibility of publish- 
ing a paper that looks disident, even if it could be correct. 

I consider that fermions, such as the electrons, satisfy the Quantum postu- 
lates because they are packets with very specific energy levels, which amounts 
to say that they have some very specific frequencies and wavelengths. For 
example, Millikan's experiment detected only discrete levels of energy, cor- 
responding to the finite number of electrons in the oil drops. 

On the other side, I consider that the photons (and in general bosons) 
really do not have discrete energy levels. They are energy quanta (in the 
sense that they are packets) that can have any positive continuous frequency. 

It is a suspicion of the present author that Planck was misled by some 
incorrect experimental results by H. Rubens, who visited Planck in October 
1900, and others. As Planck himself notes: "Whilst for small values of the 
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energy and for short waves, Wien's law was satisfactorily confirmed, note- 
worthy deviations for larger wavelengths were found, first by O. Lummer and 
E. Pringsheim, and finally by H. Rubens and F. Kurlbaum" [20] • See, for 
example, the strange curves transcribed as figure 2 in the paper by Abra- 
ham Pais and attributed to Rubens&Kurlbaum. All the curves (Wien, 
Rayleigh, Thiesen, Lummer and the observed values) are clamped around 
two fixed points at °C and 1000 °C, with a difference of up to 30% between 
the observed values and the Wien curve at around 500 °C. Besides, nobody 
would accept by now that the observed values are much better modeled by 
the Rayleigh- Jeans than by the Wien's equation, as shown in such figure. 

The "satisfactory confirmation" of Wien's law was given, evidently, by 
the careful measurements, "very good ones" [21], by Friedrich Paschen, who 
presented contributions to the Prusian Academy of Sciences between 1987 
and 1899 that showed a "perfect agreement with Wien's theory" [3]. 

As Planck felt the need to take into account those new experimental re- 
sults and "improve" Wien's law, "Planck interpolated between the laws of 
Wien and Rayleigh- Jeans but found that he could only derive a satisfactory 
equation using what he thought was merely a mathematical trick, namely 
that light is only emitted in packets. He did not for many years believe that 
these packets, known as quanta, corresponded with reality." 02] 

Nothwithstanding the fact that quantum theory, a theory generally ac- 
cepted by now, was originated in Planck's theory, it will be suggested in this 
paper that the Planck's blackbody radiation equations might be incorrect. 

2. THE PLANCK'S EQUATIONS 

The temperature of a body is a measurable effect furnished by the statis- 
tical average of the energies of the enormous amount of photons emitted by 
that body. However, this is not the practical way to compute the temperature 
because we don't know the temperature of the individual photons. 

I have two objectives in this and the following section: first, summarize 
some of the equations used in statistical thermodynamics and, second, at- 
tempt to put in evidence that the temperature, of a blackbody in particular 
but in general the temperature of any body, can profitably be interpreted in 
Physics as a certain frequency of the photons emitted by that body. 

First, the energy density radiated by a black body per unit volume, for 
a given wavelength A and absolute temperature T, is expressed by Planck's 
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radiation formula: 



P(A,T) = 



■ IT ■ h ■ c 



joule 



(1) 



where k is the Boltzmann constant k = 1.3805 -10 ■ joule/ K) 
h is Planck's constant and c is the speed of light in vacuum. 

The total energy radiated by the body at a given temperature is given by 
Stefan-Boltzmann law, which is obtained integrating the previous equation 
(1) for all values of wavelength from zero to infinity: 

E = a-T^ (2) 
where a is a constant given, in terms of more basic constants, by: 

'■^'■"^ (3) 



15 ■ /l3 . c3 

The numerical value of this constant is: a = 7.562473 ■ 10~^^ ■ ^""^^4 

The same resulting equation (2) and constant "a" are also obtained by 
integrating, for all values of the frequency, from zero to infinity, the expres- 
sion of Planck's radiation formula expressed in terms of frequency instead of 
wavelength, which is: 



■ IT ■ h - 



joule 



(4) 



exp {^) - 1 

Let us call maxnu the argument of the exponential function that appears 
in this equation: 



maxnu 



h ■ V 



(5) 



Let us replace this definition in (4) and solve for maxnu for the point at 
which equation (4) predicts the maximum radiation per unit frequency. The 
maximum is obtained when maxnu is equal to 3*(l-exp(-ma2;nn)). That is, 
when maxnu has the value 2.821439. . . 



6 



On the other hand, if we call maxlambda the argument of the exponential 
function of equation (1): 

h ■ c 

maxlambda = ^ j, ^ (6) 

Then, the point of maximum radiation per unit wavelength, computed 
from equation (1) by replacing in it this definition (6), is obtained when 
maxlambda is equal to 5*(l-exp(-maa;/amMa)). That is, when maxlambda 
has the value 4.965114... 



Let us remember the well-known relation between the wavelength and 
frequency of a given photon, which says that its product is equal to the 
speed of light: 

A • 1/ = c (7) 

if c/A from this expression is replaced by v in the right hand side of 
equation (6) then both variables, maxnu and maxlambda are proved to be 
equal to exactly the same common symbolic expression. However, the said 
common expression does not have a unique numerical value, but two values, 
given respectively by the variables maxnu and maxlambda computed above. 

If we were expecting, intuitively, that the frequency and wavelength of 
the photons emitted at the point of maximum radiation density should be 
the same, whether they are computed with equations (5) or (6), this proves 
to be wrong. 

The reason why these two constants, maxnu and maxlam,bda, are numer- 
ically different in current Physics is because there are two kinds of photons 
that participate in the computation of each one of these variables. So, the 
problem remains as to which one of them, both, or none, is the correct pho- 
ton emitted at the physical point of maximum radiation. As was said before, 
the formula (6) is usually used to compute the properties of the photons, 
but with no clear explanation as to why the different value produced with 
equation (5) is incorrect or does not correspond to the same maximum. 

In principle, either of the last two relations, (5) or (6), could be used to 
compute the values of the frequency or wavelength of the photons emitted 
when the radiation has its maximum value. Namely, solving (6) for wave- 
length: 

1 h 

Xmax = ; — — — • c • — (8) 

maxlambda ■ T k 
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Consequently, from this equation, we obtain: 



Xmax - T = b 



(9) 



This is Wien's law of displacement that allows to compute the wavelength 
of the photons emitted at the point where the density of radiation per unit 
wavelength is a maximum, by dividing the constant b by the absolute tem- 
perature of the body. 

b is called Wien's displacement constant and has units of length multiplied 
by temperature: 

b = ^"^ , , = 0.0028980841 ■ m ■ K 

k ■ maxlambda 

On the other hand, solving equation (5) for the frequency, we find an 
unnamed equation that is never used to compute the frequency of the same 
photons: 



Ok, not never but at least seldom used, because equation (8), which is 
derived from equation (6), is the most commonly used equation. However, 
as examples of the use of equation (10), which is derived from equation (5), 



The results obtained with these equations are compared in the following 



Obviously, the frequency and wavelength provided by equations (8) and 
(10), do not belong to the same kind of photons, emitted at the point of 
maximum radiation. 

3. THE TEMPERATURE IS GIVEN BY THE FREQUENCY 
OF PHOTONS 

One of the strongest hints, to accept as true that the temperature is 
frequency is obtained by replacing the wavelength from Wien's law, equation 
(9), in equation (7) and solving for T: 



Umax = maxnu ■ T ■ — 

h 



(10) 



see i, 0, dn, [H], uni, m- 



section. 



T = Umax ■ - 
c 



(11) 
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That is, the temperature of a body is equal to a frequency multiphed by 
a constant that simply converts units. Usually we know the temperature and 
wish to know the frequency, so let us define the constant kel to convert the 
temperature in Kelvin degrees to Hertz as: 

c Hz 

kel = - = 1.0344505 ■ lO" 

K 

Once the frequency is computed, multiplying this constant by the tem- 
perature of the body, we are able to solve equation (7) for wavelength. Al- 
ternatively, divide the constant h by the temperature to obtain wavelength. 

For example, to a temperature of 5530 °K (temperature of the sun, com- 
puted by Wien in his Nobel lecture |ST]), correspond photons with wavelength 
of visible light of about 524 nm (Wien reports 532 nm. Verify in jH] that this 
corresponds to a sun of green color, assuming that the photons emitted at 
the point of maximum radiation have some relation with the sun's color) and 
frequency around 5.72*10^'^ Hz. Compare this with the value of 3.25*10^^ Hz 
obtained with the use of equation (10). 

The product of wavelength and frequency of the two kinds of photons 
for which the radiation is a maximum (obtained from the two density of 
radiation equations, (8) and (10) ) is about 0.568 of the speed of light (try 
for example multiplying the values reported in ^3], JH]), which is the same 
proportion as between maxnu and maxlamhda. This proves that they are not 
the same kind of photon. 

More evident and dramatic becomes the difference if we solve equation 
(10) for the temperature that corresponds to the wavelength used by Wien, 
532 nm. This gives around 9586 °K. This means that there is a difference of 
more than 4000 °K between the sun's temperatures predicted by two suppos- 
edly valid equations of current thermodynamics. I wonder how the physicists 
justify this sort of inconsistencies or are blind to them. In the particular case 
of the sun's temperature, for example, why is it usually computed with equa- 
tion (8) instead of with equation (10)?. Several references, including the same 
Wien , insist that the point of maximum radiation of the sun is somewhere 
in the visible range, mainly between 475 and 532 nm, which corresponds to 
a sun of blue, turquoise or green color!... physicists have a lot of things to 
explain. Nevertheless, if we assume as known the sun's temperature of 5530 
°K and compute the corresponding wavelength for the maximum, again with 
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the help of equations (10) and (7), we find 922 nm, which is in the infrared, 
well outside the visible range. 

It should be noted that, in practice, neither of these equations seems to be 
used. For example in ^17j the authors report that equation AT = 3670/im ■ K 
is a "more useful displacement law than Wien's" (without explaining why it 
is more useful or how it is obtained). This equation is far away from both 
the 2898.08/xm ■ reported earlier in this paper, and the 5100/im ■ pro- 
duced respectively by equations (9) and (10), but much closer, within 2%, to 
the value that can be computed with the theory proposed in the following: 
X ■ T = h ■ c/4 ■ k = 3597.33/im ■ K. With this equation and with the more 
updated temperature for the sun of around 6000 °K we are able to compute a 
wavelength of 600 nm, which corresponds to a nice orange sun (see again [Hj). 

An interesting application of the constant kel converts Boltzmann's con- 
stant into Planck's constant. This means that Boltzmann's constant (by the 
way, never written as such by Boltzmann but by Planck, and whose currently 
accepted value might be a little greater than needed for the theory being de- 
veloped here) is redundant and superfluous if we use only frequencies instead 
of temperatures in our computations. 

4. DEVELOPMENT OF A NEW EQUATION FOR RADIA- 
TION 

The criteria used to develop a new equation for radiation, which hope- 
fully will provide an expression for the radiation law that avoids the above- 
mentioned inconsistencies, are the following: 

a. - The product of wavelength by frequency for a certain kind of photons 
should give the speed of light. 

b. - The point of the maximum radiation, expressed either in terms of 
wavelength or of frequency should be the same. 

c- The photon energy is assumed to have a continuous spectrum. This 
means that the photon energy is quantized but does not have discrete values. 

The conditions "a" and "b" are clear and intuitive and the last one was 
explained in the introduction in the sense that the physical particles show 
quantum behavior; this means that their "size" (wavelength) and energy (fre- 
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quencies) are of some fixed values. But the experimentalists have not been 
able to find discrete vahies for the corresponding photon variables. Rather, 
the values of photon energy and wavelength have a seemingly infinite and 
continuous range. 

In an open system, the radiation process emits photons to the space, 
which are not absorbed again by the system, so the mass or temperature of 

the system is reduced. 

But what we will be considering in the following is the case of a closed sys- 
tem, which is in thermal equilibrium. 

Consider an individual photon emitted by a black body. This photon 
bumps with other particles and loses and gains energy in this process, before 
it is completely absorbed or eventually gets its opportunity to get out through 
the small window of the blackbody for its energy to be measured by us. In 
fact this opportunity is so small that we are able to consider the blackbody 
as a closed system in thermodynamic equilibrium. 

Let us try to determine what the distribution of the energies of the 
individual photons is. 

The statistical weight (number of possible states of a system) of a given 
(discrete) distribution is 

ATI 

When the distribution is continuous, each partition i has only one particle 
in it and hence each factorial in the denominator is the unity. Moreover, the 
factorial in the numerator must be interpreted as the Gamma function, which 
is the continuous equivalent of the factorial. 

Consequently, the entropy S, that is just the statistical weight fl, is de- 
fined as 

S = Q = T{N) = T{p) (13) 

where the change from N to the new parameter p is to try to understand 
that now we have continuous and not integer values. 

In the following we will need to study the actual entropy distribution of these 
states among the different values of a certain variable, x, such as frequency, 
and not merely the total number of states. Therefore, the appropriate form 
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of the entropy will be: 

S{x) ^r{p,Q,x) (14) 

This provides the cumulative number of states for the variable between and 
some value x. The cumulative normalized entropy becomes: 

^,^,_ r(p,o,x) _ r(p)-r(p,x) _^ rjp, x) 

Also, to convert the units of the variable a; to a non-dimensional form, we 
include a constant. A: 

S(.) ^ 1 - (16) 

We assume that a system in a state of equilibrium is changing between 
microcanonical states in a reversible form. In such case, remember that for 
reversible processes the change of entropy is given, by the Clausius' formula, 
as proportional to the change of heat: 

dQ{x) 

dS{x) = -jr- (17) 
where T is the average temperature of the system. 

Consequently, replacing (16) in (17) and taking the derivative, we ob- 
tain the equation that describes the probability distribution of energy in a 

continuous system, as is assumed the blackbody. More specifically it also 
describes the specific heat of compounds, where the low-temperature zone 
corresponding to the so-called Schottky anomaly has this distribution: 

dQ{x)_^ dS{x)_^ y-x^-^ 
dx dx Exp{X ■ x) ■ r(p) 

The previous is the most simple derivation, which proves that the energy 
distribution per unit a; is a Gamma distribution. 

An easier and profitable understanding of the following section is achieved 
if we interpret the energy distribution as a probability distribution function 
(pdf) and the entropy as the corresponding cumulative distribution function 
(cdf ) . It should be noted that "according to Boltzmann, entropy is a measure 
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for physical probability" PB] . 

5. DERIVATION FROM FIRST PRINCIPLES 



Let us assume a continuous random variable "X" that represents some 
physical quantity such as the amount of money or energy possessed by agents 
or actors in a certain system. 

The probability that X lies in a small interval (x, x+Ax) is 

P[x<X <x + Ax]= F{x + Ax) - F{x) = E(^±^LlIi^Ax (19) 

If Aa; is small and F{x) is differentiable then the following approximation 
holds: 

P[x < X < X + Ax] = F'{x)Ax (20) 

in which the prime denotes the derivative of F{x). Whenever it exists, 
the derivative of F{x) is called the probability density function (pdf) of "X" 
and is denoted by f{x). Thus, 

/(.) ^ '-^ (21) 

The units of the pdf are inverse to those of the random variable X. 
Therefore, given the pdf, the distribution function of X is computable as 

/X 
f{u)du. (22) 
-oo 

This is the probability that X is less than x. F{x) is called the cumulative 
distribution function (cdf): 

P[-oo <X <x]= F{x) (23) 
Since the cdf is a non-decreasing function, the pdf is always nonnegative: 

f{x) > (24) 

In the following we are going to argue that, in a model in which the 
agents interact exchanging the quantity X, the equilibrium distribution of 
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this variable is a Gamma distribution. 

We define a kinetic continuous model in which, for facihty of analysis, 
consider a finite number N of agents that interact with each other by "trad- 
ing" or "colliding" between them, and exchanging random amounts of the 
quantity X, generally referred to as "money" (or energy) in the following. 

The evolution of the system is then carried out in the following way. At 
every time step two agents i and j are extracted randomly and an amount 
of money Ax is exchanged between them, 

x'j = Xj + Ax. (25) 

It can be noticed that in this way the quantity X is conserved after each 
transaction, x[ + x'- = Xi + Xj, where x\ and x'j are the agent's wealth after 
the transaction has taken place. 

The proportion of the variable x (money) owned by the system at a 
certain point x is: 

xf{x) 

Then, the amount of money, a, corresponding to a small interval Ax 
around x, is computed as the product of x by the probability density of 
occurrence of such amount at that point and by the width of the interval: 

a = xf{x)Ax (26) 
Similarly, for the point x+Ax: 

b^{x + Ax)f{x + Ax) Ax (27) 

Let us call delx the amount of money given by the agent at {x+ Ax) to 
the agent at x. 

If the system is in a state of statistical equilibrium, the amount of money 
possessed by the agent at x after the transaction, let us call a\ must be 
equal to the amount of money possessed by the (wealthier) agent at {x+ Ax) 
before the transaction, that is b, minus the amount given by this agent in 
the transaction, that is delx; and vice versa, &' is equal to the amount of 
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money possessed by the agent at x, that is a, plus the amount received from 
the other agent, that is delx. In the form of equations this means: 



a' — b — delx 
b' — a + delx 



(28) 
(29) 



But, at cquiUbrium, the amounts of money before and after the transac- 
tion must be equal, therefore a' = a and 6' = b, so the two previous equations 
are in fact the same one. 

These equations reflect the fact that in a one-dimensional, completely elastic, 
head-on, collision, classical mechanics predicts that all the energy of one of 
the agents is transferred to the other agent and vice versa. The amoimt of 
"saving" in the collisions is dependent, in general, on the number of degrees 
of freedom. 

If we assume that there is a first agent at position x that is in equilibrium 
with the second consecutive element at {x+ Ax), and if this second element 
is in equilibrium with the third, then the first is also in equilibrium with the 
third and thus, successively, we conclude that two agents at any distance will 
be in equilibrium. 

Next, the amount of money traded between these two agents is assumed 
to be composed of two terms. 

The first term is proportional to the difference between the amounts of 
money possessed by the two agents, AF{x): 



Where p is the constant of proportionality. This equality is justified 
above: the derivative of the probability distribution of a variable is equal to 
its pdf. 

Actually F{x) is probability; hence, in order to fix the units it must 
be multiplied by the amount M of total money of the system, because the 
variable interchanged in the process is assumed to be money; however, this 
amount can be seen that is simplified in the following equations. 




(30) 
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The physical interpretation of equation (30) and the objective of all this 
exercise is the following: if both agents have the same amount of money then 
the difference between these amounts is zero and there is no money inter- 
change since the system would be in equilibrium. But if we think that one 
of the agents (e.g. a photon) has more money (energy) then it can give more 
money in each transaction (energy in a collision) , in proportion to its greater 
wealth (energy difference). 

In the previous expression, the constant of proportionality, is a func- 
tion of some given positive numbers usually extracted from a uniform distri- 
bution between and 1. Since the number of transactions is very large we 
can safely assume that the net effect, if p is a random variable instead of a 
constant, is equivalent to its average. 

The second term is the amount "saved" by the wealthier agent and is 
proportional to the absolute amount of money possessed by that agent: 

Xh (31) 

Or, expanding b 

\{x + Ax)f{x + Ax)Ax (32) 

where the variable A is a function of the saving propensity of the agent 
and has inverse units to those of the variable x. 

In summary the amount of money given by the wealthier agent to the 
poorer one is 

delx = (p f{x)Ax - X{x + Ax) f{x + Ax) Ax) Ax (33) 

Replacing (26), (27) and (33) into (29) gives 

{x + Ax)f{x + Ax) Ax = xf{x)Ax + {pf{x)Ax — X{x + Ax)f{x + Ax) Ax) Ax 

(34) 

Simplifying one Ax and transposing the first term in the right hand side 
to the left 

(x + Ax)f(x + Ax) - xf(x) = pf{x)Ax - X{x + Ax)f{x + Ax) Ax 
The first term is expanded to 

xf{x + Ax) — xf{x) + Axf{x + Ax) 
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Dividing both sides by (xAx) 
fix + Ax) - fix) ^f(. + A.) ^ £^(,)_,.^(,^^,). A . A. . /(. + A.) 



Ax XX X 

Finally, taking the limit when 

Ax ^ we obtain the differential equation 

^^ = (P-l)^-^-/(^) (35) 

or 

fix) (p-i) 

w^^— — ' '''' 

whose solution is 

/(a;) =c-xP-^e-^-" (37) 

where c is some constant of integration that is computed integrating the 
distribution and normalizing to 1 and whose result is 

'''' 

The distribution thus obtained, /(x), is the powerful Gamma distribution. 

If the value of the variable p is particularized to an integer value then 
this distribution converts into the Erlang or Poisson distributions. 

If the variable p has the value 1, the Gamma distribution converts into the 
negative exponential distribution, also called Boltzmann, Gibbs, Boltzmann- 
Gibbs or simply exponential law. 

In order to have an expression without the exponential but only with the 
power term we need that the parameter A be zero. But, in that case, it can 
be proved that this distribution does not converge. 

However, in the area of econophysics it is very common the use of the 
so-called Pareto or power law distribution, although it is obvious that it is 
not a probability distribution j22I. It should be profitably replaced using the 
Gamma distribution with the proper parameters. 

If a quantity x is distributed according to the Gamma distribution then 
its average is 

(x) = ^ (39) 
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As a result, if we keep the average of x equal to the unity then p must 
be equal to A. Nevertheless the larger the value of A the smaller the variance 
since this is given by: 

The maximum of the Gamma distribution is at the position 

p — 1 1 

Xmax = ^ A ^ ' 

This means that the maximum is always towards the left of the average, 
except in the limit when A ^ oo. 

The parameters A and p used in this distribution should have in general 
values greater than or equal to 1. Accordingly, if the saving propensity of 
the agents is presupposed to be some variable between zero and 1 we will 
have to compute its inverse in order to get the needed parameter A. The 
maximum is displaced toward the left of the average precisely in proportion 
to the saving propensity. 

"The functional form of such a distribution has been conjectured to be 
a Gamma distribution on the base of an analogy with the kinetic theory of 
gases, which is consistent with the excellent fitting provided to numerical 
data" 



6. ENERGY DISTRIBUTION OF AN IDEAL GAS 

The Maxwell equation describes the energy distribution of the molecules 
of an ideal gas. This was one of the first applications of statistical methods 
in Physics. James C. Maxwell obtained it for the first time around 1857. It 
has the form 

dn _ 2ttN eI 

dE~ (^T, . kT)l Exp{§) ( 

where dn/N is the proportion of molecules with energy between E and 
E+dE. 

The Gamma distribution can be defined as the following function 
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We recover the Maxwell equation after replacing the parameters of the 
Gamma distribution with the values p=3/2 and A =l/kT: 

P(l^,E) (44) 

With these values we obtain the average energy of the particles P]: 

(E) = ^ = 3kT/2 (45) 
A 

This means that the Maxwell equation has simply been a Gamma distri- 
bution with some absolutely specific parameters {p=3/2, A =l/kT). 



In the equation for energy distribution it is possible to replace the energy 
variable everywhere by its equivalent kinetic energy, in terms of velocity 
squared. The definition suggested above provides a solid mathematical frame 
of reference for not doing that in an ad hoc fashion. The distribution of speeds 
becomes: 



And the average of v'^ is obtained from the formula for average energy, 
replacing E by the kinetic energy: 

(Im v^) = p/X = 3kT/2 (47) 

or 

{v^) = 3kT/m (48) 
7. AN APPLICATION OF THE GAMMA DISTRIBUTION 



No matter what specific values we assign to the parameters and "A" , 
they always produce a probability distribution (integral equal to unity) so 
the values selected for a given application must be given experimentally or 
by other means. 

In the previous section, for example, we used the particular values p=3/2 
and A =l/kT, with which we were able to reproduce the Maxwell equation. 
In the following we will try to find the appropriate values to describe the 
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blackbody radiation distribution. 



The first alternative is to try to obtain the exponent 3 in the variable v 
of the Gamma distribution, in order that it be equal to the corresponding 
term in the Planck distribution expressed in terms of frequencies. This gives 
the value 4 for the parameter which has the support of Wien's analysis. 

As was explained for equation (39), the value of j> fixes automatically the 
value of the variable A to the same value 4 since we can assume that the 
average frequency is the unity. 

In such case, the equivalent to the Stefan-Boltzmann constant becomes: 

4^ ■ 



Whose numerical value is: a = 5.927 ■ 10 ^ • [Wm "^K ^], given the current 
value of and is close to the current value: a = 5.668 ■ 10^^ ■ [Wm~^K~'^]. 



Normalizing the Planck distribution to unit area becomes: 

7r4(eP-^ - 1) 



(50) 



It is illustrative to compare this Planck distribution (with p=4:) with the 
Gamma distribution P(4,4,2;), in Figure 1. 

Although both distributions are close it must be noted that they are 
different. The maximum for the Gamma distribution is exactly at x=0.75 
whereas the maximum for the Planck distribution is at x=0. 70536. The 
maximum values are respectively 0.896167 and 0.875546, with around 2.3% 
percent difference between them. 

Now let us try to fit both the Planck and Gamma distributions to the 
experimental data from the cosmic background radiation, collected by the 
COBE satellite 123,0. 

The plot of the COBE data (blue dots) together with the Planck distri- 
bution with parameter p=4 and scale factors vertical=437.987 and horizon- 
tal=7. 72655 (green), and Gamma distribution with parameters p= A=4 and 
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Fig.l: Comparison between Gamma and Planck Distributions 




0.5 1 1.5 2 2.5 



Figure 1: Gamma vs. Planck 



scale factors vertical=427.909 and horizontal=7. 26667 (red), is shown in Fig- 
ure 2. The adjustment was made around the point of maximum radiation. 

Although there is a good fit of both curves there are some zones where 
we can note some discrepancy. 

Just to experiment I made the fit of a Gamma distribution with param- 
eters p=A=3.5 and scale factors vertical=45 1.258 and horizontal=7.36276 
that gives a visually almost perfect fit (Figure 3). 

The figures suggest that neither the Planck distribution nor the Gamma 
with parameters p= A=4 are perfect to adjust the COBE data. The Gamma 
distribution with p= A=3.5 is much better; however, it is necessary to ana- 
lyze more results in order to confirm or refute the possibility to change the 
parameter p to the new value 3.5, or to some other value. 

In the same sense that the Gamma distribution integrates to the unity for 
any combination of "p" and "A", both forms of the Planck distribution in- 
tegrate to the same (Stefan-Boltzmann) expression. But this does not mean 
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that both equations are vahd, as it is not true that any form of the Gamma 
distribution is adequate for expressing the energy distribution of gases, pho- 
tons, money, etc. We have to find the correct values for "p" and "A" . 

In general, "p" depends on the degrees of freedom and "A" depends on 
the kinetic energy of the system. 

8. CONCLUSIONS 

It is experimentally found that a black body produces a radiation distri- 
bution that has a given maximum density. Associated with this maximum 

there must be some kind of photons that carry the maximum amount of 
energy per unit time (or power). Those photons should have some precise 
photon frequency; however let us assume that such frequency is not specified 
at the moment. 

In particular, the current blackbody radiation theory seems to lack an 
explanation as to what the (frequency or wavelength) characteristic of this 
kind of photons is. 

This paper evidences that two different maxima, with two different kinds 
of photons, are computed, corresponding to the two Planck's equations. One 
of the equations describes the blackbody radiation energy distribution per 
unit wavelength, expressed as a function of wavelength, and the other equa- 
tion describes the blackbody radiation energy distribution per unit frequency, 
as a function of frequency. 

The variables maxnu and maxlambda, defined in section 2, can be used 
to compute the values for the two corresponding kinds of photons. 

The two different kinds of photons proposed by current physics seem to be 
just mathematical values necessary to compute the said maxima for the two 
corresponding Planck's radiation distribution densities, but it is not proved 
that either one of these values necessarily corresponds to the photons which 
carry the maximum radiation power. 

In this paper it was proved that if we assimilate the wealth of the persons 
in a population to continuous photon energy, then the probability distri- 
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bution, for a population in thermodynamic equilibrium, of the number of 
persons with increasing wealth, seems to follow a Gamma distribution. Now, 
if we multiply the number of persons with a given wealth by its correspond- 
ing wealth we obtain the amount of money held by that population group. 

Similarly, the number of photons produced by blackbody radiation seems 
to follow another Gamma distribution. If we multiply the number of pho- 
tons carrying a given energy by its corresponding energy (or frequency by 
Planck's constant) we obtain the amount of energy held by that kind of pho- 
tons (assume per unit of time). 

Consequently, if the temperature variable of equation (18) is replaced by 
the corresponding photon frequencies this should provide the actual power 
distribution. 

It should be realized that the quotient of the speed of light divided by the 
photon wavelength is the same variable as photon frequency; this means that 
there should be only one probability distribution. In other words, there is no 
way for frequency and wavelength to be interpreted as independent variables 
with different and independent probabihty distributions. 

In current Physics now and always remains the same centuries-old spirit 
that nothing can be advanced in Physics and that the old theories are cor- 
rect. Fortunately there always are some stubborn people that do not believe 
in dogmas and tries to understand and explain the behavior of Nature. For 
example, it is known that Planck himself "was attracted to Physics despite 
having been advised that nothing essentially remained to be discussed in this 
branch of learning" or that "in this field, almost everything is already 
discovered, and all that remains is to fill a few holes." [321 

Although a good fit to the Cosmic Background Radiation data was shown 
in this paper, it should be taken with skepticism. For the previous version 
of the present paper I was convinced that the COBE experiment was a good 
testbed for my equations, but later I investigated a little more and found 
that such experiment assumed the prior validity of the Planck distribution 
and only tried to detect discrepancies with respect to such distribution. As 
explains Dr. Dale Fixsen, former director of the COBE project, in a private 
communication: "The FIRAS measurement was a differential measurement. 
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The measurement was not the spectrum of the CMBR but the difference 
between the CMBR spectrum and the blackbody spectrum provided by the 
on-board external cahbrator. The display of the spectrum implicitly assumes 
the Planck spectrum is a good description of the blackbody spectrum." fT^ . 
Also, Dr. John Mather explained to me that "the CORE results are based 
on assuming that the Planck formula is correct, so you can not logically use 
the CORE results to test your theory." [TH] . 

Even though such project cannot be blamed for having used what at the 
moment could have appeared as the correct theory, I expect that new exper- 
iments be carried out to determine whether the here proposed distribution 
is a better physical model. 

The Gamma distribution is an all pervasive function that, although with- 
out any theoretical support, has been applied successfully in numerical ap- 
plications such as the distribution of money among several agents (econo- 
physics) or could profitably be applied to describe the energy and speed 
distributions in ideal and real g 3iSGS clS 8b function of temperature (or, better, 
as a function of frequency), distribution of population in cities, etc. While 
the origin of the Gamma distribution involves continuous variables its range 
of applicability should not be strictly limited to them. 

In relation to its practical use, Wien's equation is a very good one. As a 
Physics teacher says: "Wien's law does not look too bad either, we can cal- 
culate an approximation for these constants A and f3, and the new formula 
fits everywhere" ||6j. 

We have to distinguish among the new form for entropy, proposed in 
the present paper as the incomplete Gamma function, and the probabilistic 
Gamma distribution, which is obtained by the derivative of the incomplete 
Gamma function. 

Note that the Euler's Gamma function satisfies: 



^0 

This expression is the equivalent, for the Gamma distribution, of the 
partition function defined in classical thermodynamics for the Roltzmann 
distribution. 




(51) 
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